Signal detection

ABSTRACT

A transmitted signal and a received signal are combined and the combination is expected to determine whether or not a target signal is present. Either the transmitted signal or the received signal is combined with an auxiliary signal containing a range of frequencies corresponding to an anticipated Doppler shift, so that an output of the combined transmitted and received signal will be present only if a target signal exhibiting a Doppler shift within the anticipated range is present. The auxiliary signal preferably comprises finite-duration signal portions of different types so as to provide a substantially uniform frequency response throughout the selected range.

[0001] This invention relates to a method and apparatus for detecting signals. The invention is especially, but not exclusively, applicable to systems which may be subject to an unknown Doppler frequency shift affecting modulated coherent signals used for communication and ranging purposes, and particularly substantially continuous signals modulated by random or chaotic waveforms utilized by obstacle detection or collision avoidance systems.

[0002] Such systems are designed to operate in multiuser and often hostile environments, and are intended for a wide range of applications, including automotive applications, industrial robotics and unmanned vehicle navigation.

[0003]FIG. 1 is a block diagram of a typical microwave obstacle-detection system. The system comprises a signal generator 1 that produces a substantially continuous waveform x(t) with suitable bandwidth to provide required range resolution. The waveform x(t) may be deterministic (periodic or aperiodic), chaotic or purely random.

[0004] The system also has a microwave oscillator 2 that generates a sinusoidal signal with required carrier frequency, a modulator 3 that modulates one or more of the parameters (such as amplitude, phase, or frequency) of the carrier signal with the modulating waveform x(t), a power amplifier (PA) 4 that amplifies the modulated carrier signal to a required level, a microwave transmit antenna (TA) 5 that radiates an electromagnetic wave representing the modulated carrier signal towards an obstacle 6, a microwave receive antenna (RA) 7 that receives an electromagnetic wave reflected back by the obstacle 6, an input amplifier (IA) 8 that amplifies a signal provided by the receive antenna (RA) 7, and a coherent demodulator 9 that processes jointly the reference carrier signal supplied by the oscillator 2 and a signal supplied by the input amplifier (IA) 8 to reconstruct a time-delayed replica y(t) of the modulating waveform x(t).

[0005] The modulating waveform x(t) and its time-delayed replica y(t) are then processed jointly during a specified time interval by a suitable processor 10, such as correlator, to produce an estimate of the unknown time delay that is proportional to the distance (range) between the system and the obstacle 6.

[0006]FIG. 2 shows an example of a correlation function of a synchronous random binary waveform.

[0007] When there occurs a relative movement between a ranging system and an obstacle of interest, an electromagnetic wave reflected back by the obstacle and received by a coherent system will exhibit a Doppler frequency shift. The value ω_(D0) of this (angular) frequency shift can be determined from: $\omega_{D0} = {\frac{2\upsilon_{0}}{c}\omega_{0}}$

[0008] where ν₀ is the radial speed (i.e., the range rate) of a relative movement between the system and the obstacle and ω₀ is the (angular) carrier frequency of a transmitted electromagnetic wave having velocity c.

[0009] The signal reflected back by a moving obstacle can be expressed as:

z _(r)(t)=γx(t−τ ₀)cos [(ω₀+ω_(D0))(t−τ ₀)+θ]

[0010] where γ is the round-trip attenuation, τ₀ is the delay corresponding to the range D, ω_(D0) is the Doppler frequency shift, and θ is an unknown constant phase shift. Strictly speaking, the value of τ₀ cannot be constant for a nonzero Doppler frequency ω_(D0). However, in most practical cases it is assumed that the small changes in range D, of the order of the wavelength of the carrier frequency, cannot be discerned when determining a short-time time-delay estimate. Such assumption justifies decoupling delay and Doppler frequency measurements.

[0011] A baseband signal corresponding to the received signal z_(r)(t) has the form

y(t)=x(t−τ ₀)cos(ω_(D0) t+φ)

[0012] where φ is an unknown phase shift.

[0013] The correlator performs the following operation ${R_{x\quad y}(\tau)} = {\frac{1}{T_{0}}{\int_{0}^{T_{0}}{{x\left( {t - \tau} \right)}{x\left( {t - \tau_{0}} \right)}{\cos \left( {{\omega_{D0}t} + \phi} \right)}{t}}}}$

[0014] where the integral is evaluated for a plurality of hypothesized time delays τ_(min)<τ<τ_(max). When the observation interval T₀ is much shorter than one period (2π/ω_(D0)) of the Doppler frequency, the value of cos(ω_(D0)t_(m)+φ) is almost constant during T₀, so that the correlation integral can be approximated by ${R_{x\quad y}(\tau)} \approx {\frac{1}{T_{0}}{\cos \left( {{\omega_{D0}t_{m}} + \phi} \right)}{\int_{0}^{T_{0}}{{x\left( {t - \tau} \right)}{x\left( {t - \tau_{0}} \right)}{t}}}}$

[0015] where the time instant t_(m) is taken in the middle of the observation interval T₀.

[0016] When the correlation integral is calculated repeatedly for successive short processing intervals, each of duration T₀, the sequence of observed correlation functions may be represented by the plot of FIG. 3. The rate of change (in time) of the correlation function R_(xy)(τ) will correspond to the Doppler frequency ω_(D0). The value of this frequency can be determined by applying some suitable form of spectral analysis.

[0017]FIG. 4 is a block diagram of a conventional correlator, comprising variable delay line 11, multiplier 12 and integrator 13, followed by a spectrum analyser 14. When the total number of successive short processing intervals is large enough, the frequency spectrum S(ω), observed at the output of the analyser, will exhibit a pronounced peak at the Doppler frequency ω_(D0).

[0018]FIG. 5 is a block diagram of a multichannel correlator that uses a tapped delay line 15 to cover the entire interval of hypothesized time delays, τ_(min)<τ<τ_(max), in J steps of unit delay Δ using delay circuits 16. Multipliers 17 process the delayed signals wiith y(t) and provide the outputs to integrators 18. The required spectral analysis may be performed by a digital processor 19 implementing the Discrete Fourier Transform (DFT). The principle of operation is thus similar to that of FIG. 4.

[0019] The systems shown in FIGS. 4 and 5, and also other similar known systems, attempt to approximate the correlation integral with two parameters of interest, delay time τ and Doppler frequency ω_(D) by combining a time sequence of correlation functions determined over a number of relatively short observation intervals T₀. Because of this approximation, such systems can only provide suboptimal solutions to the problem of joint estimation of time delay and Doppler frequency.

[0020] Another prior art technique involves determining values of the following correlation integral ${R_{x\quad y}\left( {\tau,{\omega_{D};\phi}} \right)} = {\frac{1}{T_{0}}{\int_{0}^{T_{0}}{{x\left( {t - \tau} \right)}{\cos \left( {{\omega_{D}t} + \phi} \right)}{y(t)}{t}}}}$

[0021] for the entire interval of hypothesized time delays

τ_(min)<τ<τ_(max),

[0022] and also for the entire interval of hypothesized Doppler frequencies

ω_(Dmin)<ω_(D)<•_(Dmax);

[0023] the unknown phase φ should also be varied over the interval (0, 2π). In this case, values of the correlation integral, determined for some prescribed range of argument values (τ, ω_(D)), define a two-dimensional correlation function. The specific values of arguments τ and ω_(D), say τ₀ and ω_(D0), that maximise the two-dimensional correlation function provide estimates of unknown time delay and unknown Doppler frequency.

[0024]FIG. 6 is a block diagram of a suitable system capable of determining values of the correlation integral for some specified time delays τ and Doppler frequencies ω_(D). The system is formed by a correlator like that of FIG. 4, including variable delay line 11, multiplier 12 and integrator 13. However, the signal x(t) is first applied to a multiplier 20 which multiplies it by a variable Doppler signal cos(ω_(D)t+φ).

[0025] The procedure of adjusting the unknown value of phase φ can be avoided by suitably combining the values of the following two integrals: $\begin{matrix} {{R_{x\quad y\quad c}\left( {\tau,{\omega_{D};\phi}} \right)} = {\frac{1}{T_{0}}{\int_{0}^{T_{0}}{{y(t)}{x\left( {t - \tau} \right)}{\cos \left( {\omega_{D}t} \right)}{t}}}}} & (1) \\ {{R_{x\quad y\quad s}\left( {\tau,{\omega_{D};\phi}} \right)} = {\frac{1}{T_{0}}{\int_{0}^{T_{0}}{{y(t)}{x\left( {t - \tau} \right)}{\sin \left( {\omega_{D}t} \right)}{t}}}}} & (2) \end{matrix}$

[0026] evaluated for the entire interval of hypothesized time delays

τ_(min)<τ<τ_(max)

[0027] and the entire interval of hypothesized Doppler frequencies

ω_(Dmin)<ω_(D)<ω_(Dmax).

[0028] In many industrial and automotive applications, systems capable of detecting moving obstacles and employing substantially continuous modulated microwave (or other coherent radiation), developed in accordance with the prior art, will be too complicated and also too expensive.

[0029] Furthermore, obstacle-detection systems designed for multiuser environments will preferably exploit some suitable form of random or chaotic modulation which may preclude the use of conventional Doppler signal processors constructed in accordance with the prior art.

[0030] Aspects of the present invention are set out in the accompanying claims.

[0031] In accordance with a further aspect of the invention, representations of a transmitted signal and a received signal are processed together in such a manner as to determine whether or not a target signal is present (this target signal possibly representing an object reflecting the transmitted signal). An auxiliary representation, for example formed of one or more auxiliary signals, is introduced into the process. The auxiliary representation includes frequencies distributed throughout a range which corresponds to anticipated frequency modifications (for example Doppler shifts) in the target signal, as a result of which the system has response characteristics such that a target signal exhibiting frequency modifications within this range will give rise to a significant system output, whereas frequency modifications outside this range will not.

[0032] Such an arrangement avoids the need for determining the value of the frequency modification (e.g. Doppler shift) but enables the rapid detection of any target signal incorporating a frequency modification of interest, using a simple structure. The above-mentioned arrangement is thus particularly useful in applications where it is sufficient to discriminate between stationary obstacles and obstacles moving with specified velocities, and for applications in which obstacles are placed into a number of classes according to obstacle velocity, for example zero/low, moderate/average, high/very high.

[0033] Preferred aspects of the present invention are directed to the nature of the auxiliary representation and give rise to a number of advantages as explained below.

[0034] Preferably, the auxiliary representation is in the form of finite-duration signal portions comprising multiple frequencies within the range of interest, and preferably distributed throughout this range. The duration of the signal portions is selected so as to maintain orthogonality of the components of the respective frequencies. By using finite-duration portions and selecting the shapes of these portions appropriately, it is possible to extend the system response characteristics such that a significant output is produced for not only target signal frequencies equal to those of the signal components used to construct the auxiliary representation, but also adjacent frequencies.

[0035] Preferably, each signal portion is created using a window function which results in a significant frequency response for all frequencies within the range of interest. However, it is difficult to arrange for a signal portion to produce a flat response throughout the frequency range of interest.

[0036] For this reason, in accordance with a preferred embodiment of the invention, the auxiliary representation comprises different finite-duration signal portions, each comprising components having a respective set of frequencies, and arranged such that they give rise to complementary frequency response characteristics. Accordingly, by combining the results achieved using the respective different auxiliary signal portions, it is possible to get an overall frequency response characteristic which is significantly flatter than if only one type of signal portion were used.

[0037] Preferably, each signal portion comprises components of frequencies interleaved with those of the other signal portion or portions.

[0038] An indicated by FIG. 3, at any given instant the correlation function may have a low level depending upon the phase of the Doppler shift. In order to ensure a significant output irrespective of the unknown phase, the auxiliary representation preferably comprises a first representation and a second representation, the components of respective frequencies in the first representation having a quadrature phase relationship with the corresponding frequency components in the second representation. By combining the responses achieved as a result of theses two auxiliary signals, the output is substantially unaffected by the phase of the Doppler shift.

[0039] Accordingly, in the preferred embodiment, there are at least four types of auxiliary signal portions: signal portions which differ in that they comprise frequencies interleaved with frequencies of other portions, and signal portions which differ in that their frequency components are in quadrature relationship with those of other portions.

[0040] The preferred embodiment uses four types of signal portions, but it would be possible to produce a system which uses more types, while maintaining the advantages indicated above. For example, the frequency components may be split into three or more sub-sets.

[0041] It would be desirable to synthesise digitally the auxiliary representation. In order to facilitate this, the auxiliary representation preferably exhibits small values of the peak factor, i.e., the peak value divided by its root-mean square value. Minimising this enables the waveform to be represented digitally with a minimum number of bits for a specified reconstruction error.

[0042] A small peak factor can be achieved by appropriate selection of the initial phase values of the signal components for each signal portion.

[0043] In the preferred embodiment described below, the transmitted and received representations are signals which are combined by correlation. However, other arrangements are possible. For example, the received signal may be modified (e.g. multiplied) by the auxiliary representation and then delivered to a filter matched to the characteristics of the transmitted signal. This use of matched filters, which utilise a representation of the transmitted signal rather than the signal itself, is, per se, known in the art.

[0044] In another alternative embodiment, the correlator is replaced by a suitably modified version of the time-delay discriminator described in WO-A-00/39643, the contents of which are incorporated herein.

[0045] Although the present invention is being described primarily in the context of detection of obstacles which produce a Doppler shift in a reflected signal, the invention is applicable also in other areas. For example, the invention is applicable to communication systems in which the local oscillators in the transmitter and receiver may not be accurately matched; a frequency offset between the oscillators would give rise to a frequency shift similar to a Doppler shift. The invention thus permits use of equipment which has wider tolerances, and is thus less expensive, than prior art arrangements.

[0046] The invention relates both to a method of signal processing and to apparatus operating according to such a method.

[0047] Arrangements embodying the invention will now be described by way of example with reference to the accompanying drawings.

[0048]FIG. 1 is a block diagram of a typical microwave obstacle-detection system.

[0049]FIG. 2 is a correlation function of a synchronous random binary waveform.

[0050]FIG. 3 shows the effects of different values of phase of the Doppler shift in a target signal on the correlation function of a synchronous random binary waveform.

[0051]FIG. 4 is a block diagram of a conventional correlator followed by a spectrum analyser.

[0052]FIG. 5A block diagram of a multichannel correlator followed by a digital processor implementing the Discrete Fourier Transform.

[0053]FIG. 6 is a block diagram of a system for determining values of the correlation integral for various time delays and various Doppler frequencies.

[0054]FIG. 7 is a representation of the time-frequency plane divided into range-Doppler resolution cells.

[0055]FIG. 8 is a representation of the time-frequency plane divided into range resolution cells for a specified interval of Doppler frequencies.

[0056]FIG. 9 is a block diagram of a two-channel Doppler-tolerant signal processor in accordance with the present invention.

[0057]FIG. 10 is a block diagram of a multichannel correlator constructed in accordance with the prior art.

[0058]FIG. 11 is an example of an auxiliary representation in the form of signal portions for use in the processor of FIG. 9.

[0059]FIG. 12 schematically illustrates an example of operations performed by a channel combiner of the processor of FIG. 9.

[0060]FIG. 13 shows an alternative configuration of a two-channel Doppler-tolerant signal processor in accordance with the present invention.

[0061]FIG. 14 is a block diagram of a four-channel Doppler-tolerant signal processor in accordance with the present invention.

[0062]FIG. 15 is a block diagram of a quadrature Doppler generator for use in a Doppler-tolerant signal processor in accordance with the present invention.

[0063] FIGS. 16(a) and (b) and FIGS. 17(a) and (b) show examples of synthesized signal portions used to construct the auxiliary representation.

[0064]FIG. 18 shows the partial frequency response of a two-channel Doppler-tolerant signal processor employing the signal portions of FIGS. 16(a) and (b).

[0065]FIG. 19 shows the complementary frequency response of a two-channel Doppler-tolerant signal processor employing the signal portions of FIGS. 17(a) and (b).

[0066]FIG. 20 shows the combined frequency response of a two-channel Doppler-tolerant signal processor employing all four signal portions.

[0067]FIG. 21 shows ripples of the combined frequency response of FIG. 20.

[0068]FIG. 22 shows the combined frequency response of a two-channel Doppler-tolerant signal processor employing a squared cosine window.

[0069]FIG. 23 shows ripples of the combined frequency response of FIG. 22.

[0070]FIG. 24 shows a modified Kaiser window (solid line) and a squared cosine window (broken line).

[0071]FIG. 25 shows two sets, A and B, of auxiliary signals employed by two two-channel Doppler-tolerant signal processors operating in a tandem configuration.

[0072]FIG. 7 shows a representation of the time-frequency plane divided into range-Doppler resolution cells and a hypothetical obstacle occupying a single cell. Prior art arrangements are designed to establish the coordinates τ₀,ω_(D0) of the cell.

[0073] When obstacles are to be classified into only two alternative classes: obstacles moving with specified velocities and all other obstacles, the time-frequency plane will be divided into a number of resolution cells equal to the number of required range cells. An obstacle detected in any of those resolution cells will have the range specified by that cell, but its velocity may be one of the plurality of admissible velocities, corresponding to an entire interval of Doppler frequencies ω_(D), say, ω_(Dmin)<ω_(D)<ω_(Dmax).

[0074]FIG. 8 is a representation of the time-frequency plane divided into a number of range resolution cells for a specified interval of Doppler frequencies. An obstacle located at the range cell τ₀ may have a radial velocity within an interval of velocities corresponding to the entire Doppler interval ω_(Dmin)<ω_(D)<ω_(Dmax).

[0075] In order to detect obstacles moving with any of the specified velocities, the two correlation integrals, given by (1) and (2), can be replaced by $\begin{matrix} {{R_{x\quad y\quad c}\left( {\tau,{\omega_{D};\phi}} \right)} = {\frac{1}{T_{0}}{\int_{0}^{T_{0}}{{y(t)}{x\left( {t - \tau} \right)}{d_{C}(t)}{t}}}}} & (3) \\ {{R_{x\quad y\quad s}\left( {\tau,{\omega_{D};\phi}} \right)} = {\frac{1}{T_{0}}{\int_{0}^{T_{0}}{{y(t)}{x\left( {t - \tau} \right)}{d_{S}(t)}{t}}}}} & (4) \end{matrix}$

[0076] where d_(C)(t) and d_(S)(t) are especially designed auxiliary signals comprising a plurality of suitable frequency components.

[0077] The auxiliary signals dc(t) and d₅(t) should be designed in such a way as:

[0078] to contain all, mutually orthogonal within the interval To, frequency components covering the entire interval of Doppler frequencies of interest;

[0079] to contain two quadrature components at each frequency of interest, thus providing independence of unknown phase shift;

[0080] to exhibit small values of the peak factor (i.e., the peak value of a waveform divided by its root-mean-square value) required to represent an analogue waveform digitally with a minimum number of bits for a specified value of a reconstruction error;

[0081] to provide a signal processor with a flat frequency response over the entire interval of frequencies of interest, and also a required attenuation at all other Doppler frequencies.

[0082] According to a preferred embodiment of the present invention, the two auxiliary signals, d_(C)(t) and d_(S)(t), are obtained by suitably combining four finite-duration chirps, d_(CE)(t), d_(SE)(t), d_(CO)(t) and d_(SO)(t), defined by: $\begin{matrix} {{d_{C\quad E}(t)} = {{w_{E}(t)}{\sum\limits_{k = 1}^{K}{\cos \left\lbrack {{2{\pi \left( {{2k} + u} \right)}f_{p}t} + \alpha_{k}} \right\rbrack}}}} & (5) \\ {{d_{S\quad E}(t)} = {{- {w_{E}(t)}}{\sum\limits_{k = 1}^{K}{\sin \left\lbrack {{2{\pi \left( {{2k} + u} \right)}f_{p}t} + \alpha_{k}} \right\rbrack}}}} & (6) \\ {{d_{C\quad O}(t)} = {{w_{O}(t)}{\sum\limits_{k = 1}^{K}{\cos \left\lbrack {{2{\pi \left( {{2k} + 1 + u} \right)}f_{p}t} + \beta_{k}} \right\rbrack}}}} & (7) \\ {{d_{S\quad O}(t)} = {{- {w_{O}(t)}}{\sum\limits_{k = 1}^{K}{\sin \left\lbrack {{2{\pi \left( {{2k} + 1 + u} \right)}f_{p}t} + \beta_{k}} \right\rbrack}}}} & (8) \end{matrix}$

[0083] where W_(E)(t) and w_(O)(t) are suitable window (or tapering) functions, K is the number of even or odd frequency components, u is the integer frequency shift, ƒ_(p) is the fundamental frequency, and α_(k) and β_(k) are phase angles so selected as to obtain the minimum, or at least low, value of the peak factor.

[0084] The chirp d_(SE)(t) is a quadrature copy of the corresponding chirp d_(CE)(t); similarly, the chirp d_(SO)(t) is a quadrature copy of the corresponding chirp d_(CO)(t). All harmonics of the fundamental frequency ƒ_(p) are shifted (by u) and split into two sets, each comprising either even or odd harmonics. Those sets are then employed to synthesize two orthogonal groups of chirps: {d_(CE)(t), d_(SE)(t)}, and {d_(CO)(t), d_(SO)(t)}.

[0085] Window functions W_(E)(t) and w_(O)(t) belong to a large class of known and widely used window functions, such as Blackman windows, Kaiser-Bessel windows or Dolph-Chebyshev windows. The duration T_(w) of a selected window is so chosen as to exploit the orthogonality of the frequency components, and the window type and shape are mainly determined by a required frequency resolution at each Doppler frequency of interest.

[0086] In order to simplify the optimisation procedure and also to obtain symmetrical chirps, the values of phases α_(k) and β_(k) can be restricted to only two values, 0 and π. In such a case, the four chirps, d_(CE)(t), d_(SE)(t), d_(CO)(t) and d_(SO)(t) will be synthesized as follows: $\begin{matrix} {{d_{C\quad E}(t)} = {{w_{E}(t)}{\sum\limits_{k = 1}^{K}{a_{k}{\cos \left\lbrack {2{\pi \left( {{2k} + u} \right)}f_{p}t} \right\rbrack}}}}} & (9) \\ {{d_{S\quad E}(t)} = {{- {w_{E}(t)}}{\sum\limits_{k = 1}^{K}{a_{k}{\sin \left\lbrack {2{\pi \left( {{2k} + u} \right)}f_{p}t} \right\rbrack}}}}} & (10) \\ {{d_{C\quad O}(t)} = {{w_{O}(t)}{\sum\limits_{k = 1}^{K}{b_{k}{\cos \left\lbrack {2{\pi \left( {{2k} + 1 + u} \right)}f_{p}t} \right\rbrack}}}}} & (11) \\ {{d_{S\quad O}(t)} = {{- {w_{O}(t)}}{\sum\limits_{k = 1}^{K}{b_{k}{\sin \left\lbrack {2{\pi \left( {{2k} + 1 + u} \right)}f_{p}t} \right\rbrack}}}}} & (12) \end{matrix}$

[0087] where the coefficients α_(k) and β_(k) may assume only two values +1 or −1, and are so selected as to minimize the value of the peak factor.

[0088] A signal processor determining and combining the values of the two correlation integrals (3) and (4) will be referred to as a “Doppler-tolerant” signal processor. The resulting wideband Doppler compensation will make all obstacles of interest detectable (‘visible’) at the output of a Doppler-tolerant signal processor, while all other obstacles will be strongly attenuated.

[0089]FIG. 9 is a block diagram of a two-channel Doppler-tolerant signal processor, incorporating a quadrature Doppler generator QDG, according to the present invention. The processor utilizes two identical signal processing blocks, CRI and CRQ, either of which may be a conventional correlator or another appropriate signal processor. In the following, it will be assumed that each of the blocks CRI and CRQ is a conventional multichannel correlator.

[0090]FIG. 10 is a block diagram of a multichannel correlator, constructed in accordance with the prior art (see FIG. 5), which is an example of a suitable structure for each of blocks CRI and CRQ. The correlator provides J discrete values of the correlation function determined at the multiples of unit delay A.

[0091] In the two-channel Doppler-tolerant signal processor of FIG. 9, the function of Doppler compensation is implemented jointly by a quadrature Doppler generator QDG, two analogue multipliers, AMI and AMQ, and a channel combiner IQC.

[0092] A complete set of operations is performed by the two-channel signal processor during a basic processing time interval 2T_(w) split into two equal subintervals T_(w), referred to as an even cycle and an odd cycle.

[0093] A baseband signal y(t), supplied by the coherent demodulator 9 in FIG. 1 (shown here as comprising a mixer AM for mixing the input signal with a local oscillator signal, followed by a low pass filter LPF), is applied simultaneously to inputs I1 and Q1 of the multipliers AMI and AMQ, respectively. The other two inputs, I2 and Q2, of the multipliers are driven by auxiliary signals produced by the quadrature Doppler generator QDG.

[0094] During an even cycle, the generator QDG produces simultaneously two chirps: d_(CE)(t) at output I2, and d_(SE)(t) at output Q2. During an odd cycle, another two chirps, d_(CO)(t) and d_(SO)(t), are produced simultaneously at outputs I2 and Q2, respectively. When the processing interval 2T_(w) is replicated, the alternate cycles contain chirps d_(CE)(t) and d_(CO)(t) which will thus form a periodic auxiliary signal d_(C)(t). Similarly, a periodic auxiliary signal d_(S)(t) will be obtained by periodic extension of combined chirps d_(SE)(t) and d_(SO)(t). Each of these two resulting quadrature signals, d_(C)(t) and d_(S)(t), contains the same frequency components covering the entire hypothesized Doppler frequency range of interest.

[0095] The generator QDG also produces a pulse train EOC to indicate the time instants of cycle transitions, but not necessarily a particular even or odd cycle.

[0096]FIG. 11 shows an example of two auxiliary quadrature signals, d_(C)(t) and d_(S)(t), and a pulse train EOC produced by the quadrature Doppler generator QDG.

[0097] Each correlator, CRI or CRQ, processes jointly a reference signal x(t) and a signal supplied by the corresponding multiplier, AMI or AMQ, during even and odd cycles. These cycles, each of duration T_(w), are initiated by consecutive pulses of the pulse train EOC. Results presented at the output of the correlator CRI, after completing an even cycle and an odd cycle, are denoted by CE and CO, respectively. Similarly, SE and SO denote results obtained at the output of the correlator CRQ after completing an even cycle and an odd cycle, respectively.

[0098] During an even cycle, the channel combiner IQC determines a representation of

RE={square root}{square root over (CE²+SE²)}

[0099] and during an odd cycle, IQC determines a representation of

RO={square root}{square root over (CO²+SO²)}

[0100] The above representations can be determined with the use of a conventional ‘Pythagoras processor’ or another suitable processor constructed according to the prior art.

[0101] In order to simplify the structure of the channel combiner IQC, the values of RE and RO can also be determined from some suitable approximations. For example,

RE≈max{|CE|,|SE|}+0.5 min{|CE|,|SE|}

RO≈max{|CO|,|SO|}+0.5 min{|CO|,|SO|}

[0102] Other approximations can also be constructed as known in the art.

[0103] Representations of RE and RO, determined for a prescribed number L of the processing intervals 2T_(w), are suitably combined together, e.g., added or averaged, and then the resulting combined representation REO is provided at the output of the combiner IQC. For example, for a single observation interval, L=1 and REO=RE+RO. All operations performed by the channel combiner can be implemented either by a special-purpose digital processor or by a standard programmable digital signal processor.

[0104] For obstacle detection purposes, the combined representation REO is to be compared to a predetermined decision threshold. Obstacle detection will be declared in those delay (range) cells, where the threshold has been exceeded by the corresponding value of the combined representation REO.

[0105]FIG. 12 is a schematic illustration of an example of operations performed by the channel combiner IQC and the resulting detection decisions, which employ a decision threshold DT.

[0106]FIG. 13 shows an alternative configuration of a two-channel Doppler-tolerant signal processor with Doppler compensation. In this case, it is the reference signal x(t), and not a demodulated signal y(t), that is multiplied by the chirps: d_(CE)(t), d_(SE)(t), d_(CO)(t) and d_(SO)(t) produced by the quadrature Doppler generator QDG. Except for this modification, all operations and functions performed by the processor remain identical to those performed by the configuration shown in FIG. 9.

[0107]FIG. 14 is a block diagram of yet another configuration of a Doppler- tolerant signal processor utilizing Doppler compensation. The system comprises four channels: a pair of in-phase channels and a pair of quadrature channels, each with their own correlator CRIE, CRQE, CRIO or CRQO. In this case, the chirps: d_(CE)(t), d_(SE)(t), d_(CO)(t) and d_(SO)(t) are supplied simultaneously during a single cycle T_(w). Consequently, the quadrature Doppler generator has now four signal outputs, in addition to an output for a processing interval signal PI indicative of the interval over which correlation should take place. Using additional separate channels for even-cycle chirps and odd-cycle chirps reduces by half the time required for determining the combined representation REO by a suitably modified channel combiner IQC.

[0108] It is also possible to devise another configuration of the quad-channel Doppler-tolerant signal processor processor shown in FIG. 14, by exploiting the modification presented in FIG. 13.

[0109] In some applications, it may be convenient to employ a single-channel Doppler tolerant signal processor, in which the four chirps, d_(CE)(t), d_(SE)(t), d_(CO)(t) and d_(SO)(t) will be produced and used sequentially.

[0110] All the above configurations of a Doppler-tolerant signal processor, and also other similar structures, will utilize an appropriate form of a quadrature Doppler generator QDG to provide suitable auxiliary signals for Doppler frequency compensation in accordance with the present invention.

[0111]FIG. 15 is a block diagram of a suitable quadrature Doppler generator QDG. The generator comprises an N-bit binary counter BCT, driven by a clock generator CLK, two memories, CVM and SVM, two buffer registers, CBR and SBR, two digital-to-analogue converters, CDA and SDA, and a trigger circuit TRG.

[0112] The operations performed by the quadrature Doppler generator can be summarised as follows. The most significant bit (MSB) of the output of the free-running counter BCT identifies even cycles, when MSB=0, and odd cycles, when MSB=1. The MSB and remaining (N−1) bits are used jointly to define addresses of 2^(N) cells of each memory, CVM and SVM. The memory CVM contains digital amplitude samples of a chirp d_(CO)(t) at addresses with MSB=0, whereas samples of a chirp d_(CO)(t) are stored at addresses with MSB=1. Similarly, the memory SVM contains digital samples of a chirp d_(SE)(t) at addresses with MSB=0, whereas samples of a chirp d_(SO)(t) are stored at addresses with MSB=1. Consequently, each of the four chirps is represented by 2^(N−1) binary words, each comprising M bits, where the number M is chosen to attain a specified accuracy of the resulting digital representation of chirp samples.

[0113] Successive digital words appearing at the outputs of memories CVM and SVM are transferred to the corresponding buffer registers, CBR and SBR, at the time instants determined by a suitably delayed pulse train available at output 2 of the clock generator CLK. Digital words held in the buffer registers CBR and SBR are then converted respectively into analogue samples by two corresponding digital-to-analogue converters, CDA and SDA. As a result, the converter CDA produces repeatedly: 2^(N−1) analogue samples of chirp d_(CE)(t) followed by 2^(N−1) analogue samples of chirp d_(CO)(t). Similarly, the converter SDA produces repeatedly: 2^(N−1) analogue samples of chirp d_(SE)(t) followed by 2^(N−1) analogue samples of chirp d_(SO)(t).

[0114] A trigger circuit, TRG, generates a short pulse at the time instants coinciding with each transition of a binary waveform representing values of the most significant bit (MSB) of successive states of the counter BCT. A sequence of such short pulses, denoted by EOC, is used to indicate the beginning of each (even or odd) cycle. The pulse train EOC is utilized in the Doppler-tolerant signal processor, shown in FIG. 9 and FIG. 13, to define integration intervals T_(w) of the correlators, CRI and CRQ, and also to initiate various operations performed by the channel combiner IQC.

[0115] For a quad-channel Doppler-tolerant signal processor, the quadrature Doppler generator will be modified as follows. The number of bits of the binary counter BCT will be reduced to (N−1), and successive counter states will determine addresses of four memories containing digital samples of corresponding four chirps, d_(CE)(t), d_(SE)(t), d_(CO)(t) and d_(SO)(t). The output of each of the four memories will then drive a separate digital-to-analogue converter via a suitable buffer register. As a result, analogue representations of the four chirps will be available simultaneously at four outputs of such modified quadrature Doppler generator.

[0116] The modifications of the quadrature Doppler generator required for a single-channel operation will be obvious to those skilled in the art.

[0117] In order to illustrate some basic aspects of the present invention, and also to demonstrate its potential benefits, consider the following example.

[0118] Let the number of frequency components K=31 and the integer frequency shift u=4; hence, the minimum frequency and the maximum frequency are equal to 6ƒ_(p) and 67ƒ_(p), respectively. Assume that each of the four chirps is to be represented by 2048 discrete samples and that w_(E)(t)=w_(O)(t)=w(t); assume also that α_(k)=b_(k) for k=1,2, . . . 31.

[0119] It has been found that the peak factor achieves its minimum value, when

[0120] α₁=+1 α₂=+1 α₃=+1 α₄=+1 α₅=+1 α₆=+1 α₇=+1

[0121] α₈=−1 α₉=−1 α₁₀=−1 α₁₁=+1 α₁₂=+1 α₁₃=+1

[0122] α₁₄=−1 α₁₅=−1 α₁₆=−1 α₁₇=+1 α₁₈=+1 α₁₉=−1

[0123] α₂₀=+1 α₂₁=+1 α₂₂=−1 α₂₃=−1 α₂₄=+1 α₂₅=−1

[0124] α₂₆=−1 α₂₇=+1 α₂₈=−1 α₁₉=+1 α₃₀=−1 α₃₁=+1.

[0125] Application of the above coefficients has resulted in the reduction of the peak factor value from 31/{square root}(31·0.5)≈7.9 (when all components are simply added with unit coefficients) to only 1.9. It is desirable to ensure the peak factor value is less than 2.

[0126] To achieve a flat frequency response over the selected frequency range, and provide a required attenuation at all other frequencies, the discrete window function w(j) is a modified Kaiser window of the form ${w(j)} = \frac{{I_{0}\left( {\xi \sqrt{1 - m^{2}}} \right)} - 1}{{I_{0}(\xi)} - 1}$

[0127] where m=(j−1024)/1024 for j=0,1, . . . ,2047, I_(O)(·) is the modified Bessel function, and the shape parameter ξ is equal to 6.24.

[0128] Samples of the four discrete chirps, d_(CE)(j), d_(SE)(j), d_(CO)(j) and d_(SO)(j), are determined from ${d_{C\quad E}(j)} = {2{w(j)}{\sum\limits_{k = 1}^{31}{a_{k}{\cos \left\lbrack {2{\pi \left( {{2k} + 4} \right)}{j/1024}} \right\rbrack}}}}$ ${d_{S\quad E}(j)} = {{- 2}{w(j)}{\sum\limits_{k = 1}^{31}{a_{k}{\sin \left\lbrack {2{\pi \left( {{2k} + 4} \right)}{j/1024}} \right\rbrack}}}}$ ${d_{C\quad O}(j)} = {2{w(j)}{\sum\limits_{k = 1}^{31}{a_{k}{\cos \left\lbrack {2{\pi \left( {{2k} + 5} \right)}{j/1024}} \right\rbrack}}}}$ ${d_{S\quad O}(j)} = {{- 2}{w(j)}{\sum\limits_{k = 1}^{31}{a_{k}{\sin \left\lbrack {2{\pi \left( {{2k} + 5} \right)}{j/1024}} \right\rbrack}}}}$

[0129] for j=0,1, . . . , 2047.

[0130] The chirps d_(CE)(j), d_(SE)(j), d_(CO)(j) and d_(SO)(j), synthesized above, can be produced by a quadrature Doppler generator QDG, with parameters N=12 and M=12.

[0131] The synthesized chirps d_(CE)(j), d_(SE)(j), d_(CO)(j) and d_(SO)(j) are shown in FIG. 16 and FIG. 17. As seen, the chirps exhibit various forms of symmetry which can be exploited to reduce the memory of the quadrature Doppler generator QDG, at the expense of more complicated timing and control functions.

[0132] The chirps d_(CE)(j) and d_(SE)(j) are used for correlation during an even cycle to produce the results CE and SE, respectively. Those results are then combined by the channel combiner IQC as follows

RE={square root}{square root over (CE²+SE²)}

[0133] The resulting partial frequency response of a two-channel Doppler-tolerant signal processor is shown in FIG. 18.

[0134] Similarly, the chirps d_(CO)(j) and d_(SO)(j) are used for correlation during an odd cycle to produce the results CO and SO, respectively. Those results are combined by the channel combiner IQC as follows

RO={square root}{square root over (CO²+SO²)}

[0135] The resulting complementary frequency response of a two-channel Doppler-tolerant signal processor is shown in FIG. 19.

[0136] The two frequency response characteristics are superimposed (added) to produce the combined frequency response, shown in FIG. 20. As seen, the relative attenuation at out-of-band frequencies can exceed 50 dB.

[0137] The plot of the combined frequency response in FIG. 21 reveals that the ripples of the response do not exceed ±0.25 dB at all frequencies of interest.

[0138]FIG. 22 shows the combined frequency response when the modified Kaiser window is replaced a simple cosine window of the form ${w(j)} = {\cos^{2}\left( {\frac{j - 1024}{1024}\frac{\pi}{2}} \right)}$

[0139] for j=0, 1, . . . , 2047.

[0140] The cosine window provides greater attenuation of stationary obstacles (when Doppler frequency is zero); however, the frequency response is not as flat as that corresponding to the modified Kaiser window. The shape and location of ripples are shown in FIG. 23.

[0141] In all the above configurations of a Doppler-tolerant signal processor, the received signal y(t) is multiplied by chirps with nonuniform envelopes determined by a selected window. Because of this loss of signal energy, the detection performance of the system can be significantly reduced. The examination of various windows utilized in a Doppler-tolerant signal processor shows that their shape is close to that of a squared cosine function. For example, FIG. 24 shows that a modified Kaiser window and a squared cosine window almost overlap. This observation and the trigonometric identity sin² φ+cos² φ=1 suggest that the loss due to the use of windows can be almost eliminated when a received signal is processed simultaneously by two Doppler-tolerant signal processors operating in tandem configuration. While either processor employs the same set of chirps, all chirps of a duplicated set are shifted in time by T_(w)/2, where T_(w) is the window duration.

[0142]FIG. 25 shows two sets of chirps employed by a tandem arrangement comprising two two-channel Doppler-tolerant signal processors, each arranged in accordance with FIG. 9.

[0143] A more efficient system will utilize a tandem configuration comprising two quad-channel Doppler-tolerant signal processors, each developed in accordance with FIG. 14.

[0144] This technique, whereby auxiliary signals are used to investigate a primary signal, and wherein each auxiliary signal comprises finite-duration portions overlapping portions of the other auxiliary signal, is considered independently inventive. It can be applied to other areas of signal processing beyond the specific areas disclosed herein, including but not limited to arrangements in which the auxiliary signals contain components of different frequencies. It can also be extended to cover the use of more than two trains of auxiliary signals.

[0145] The correlators of the embodiments described above could be replaced by other arrangements such as a suitably modified time-delay discriminator as described in WO-A-00/39643. 

1. A method of detecting a target signal, the method involving combining a primary representation of a transmitted signal and a primary representation of a received signal, after modifying one of the primary representations using an auxiliary representation, and determining from the result of the combination whether the target signal is present, characterised in that the auxiliary representation comprises a multi-frequency signal having multiple concurrent components of different frequencies within a frequency range corresponding to possible Doppler shifts in the target signal, whereby a received signal exhibiting a Doppler shift corresponding to a said frequency produces a significant response in the combination result.
 2. A method as claimed in claim 1, wherein the auxiliary representation components have frequencies spread throughout the predetermined range of possible Doppler shifts, whereby a received signal which exhibits any Doppler shift within the predetermined range results in a significant response in the combination result.
 3. A method as claimed in claim 1 or claim 2, wherein the different frequencies are respective harmonics of a common basic frequency.
 4. A method as claimed in any preceding claim, wherein the auxiliary representation has a finite duration, the primary representations being combined throughout that duration to produce said result.
 5. A method as claimed in claim 4, wherein the duration is substantially equal to an integral number of periods of the lowest frequency component.
 6. A method as claimed in claim 4 or claim 5, wherein the auxiliary representation is formed using a predetermined window function.
 7. A method as claimed in any preceding claim, including the step of modifying one of said primary representations with at least one further auxiliary representation, the respective auxiliary representations having components of different frequencies.
 8. A method as claimed in claim 7, wherein the frequencies of the components of a first of the auxiliary representations are interleaved with the frequencies of the components of another auxiliary representation.
 9. A method as claimed in claim 7 or 8, wherein the respective results produced when the auxiliary representations are modified with a primary representation are combined.
 10. A method as claimed in claim 9, wherein the auxiliary representations give rise to complementary frequency response characteristics, whereby the frequency response of the combined results is substantially flat throughout said frequency range.
 11. A method as claimed in any preceding claim, including the step of using first and second auxiliary representations, the components of the first auxiliary representation being in quadrature relationship with corresponding components in the other auxiliary representation, the method including the step of simultaneously modifying the first and second auxiliary representations with respective versions of one of the two primary representations and combining the respective modified representations with the other primary representation to produce a result indicative of the presence of the target signal.
 12. A method as claimed in any preceding claim, including the step of providing a first set of auxiliary representations of finite duration and a second set of auxiliary representations of finite duration, and interleaving the auxiliary representations of the first set with the auxiliary representations of the second set.
 13. A method as claimed in any preceding claim, wherein the primary representations are combined by correlation.
 14. A method as claimed in any one of claims 1 to 12, wherein one of the primary representations is a signal which is combined with the other primary representation by being applied to a filter matched to said other primary representation.
 15. A method as claimed in any preceding claim, when used for detecting the presence of an object which produces said target signal by reflection of said transmitted signal.
 16. A method as claimed in claim 15, including determining from the result of said combination the range of the object.
 17. A method as claimed in any preceding claim, including the step of digitally synthesising the auxiliary representation using stored sample values.
 18. A method as claimed in claim 17 when directly or indirectly dependent on claim 4, wherein the components of the auxiliary representation have phase values at the beginning of the representation which are such as to obtain a peak factor for the representation of less than two.
 19. A method as claimed in claim 18, wherein the phase values of the components at the beginning of the representation are selected from a first value and a second value, the first and second values differing by π.
 20. A method of detecting a target signal, the method involving combining a primary representation of a transmitted signal and a primary representation of a received signal, after modifying one of the primary representations using an auxiliary representation, and determining from the result of the combination whether the target signal is present, characterised in that the auxiliary representation comprises a multi-frequency signal having multiple concurrent components of different frequencies within a frequency range corresponding to possible frequency discrepancies between the transmitted and received signals, whereby a received signal exhibiting a discrepancy corresponding to a said frequency produces a significant response in the combination result.
 21. Apparatus for detecting a target signal, the apparatus being arranged to operate in accordance with a method of any one of the preceding claims.
 22. A signal processing method in which a primary signal is modified by, at least, first and second auxiliary signals in order to investigate the primary signal, each auxiliary signal comprising successive finite-duration shaped portions having relatively low-amplitude leading and trailing parts, the portions being interleaved with, and overlapping, signal portions of the other auxiliary signal.
 23. A method as claimed in claim 22, wherein each shaped portion has a shape substantially similar to that of a squared cosine.
 24. A method as claimed in claim 23, wherein each shaped portion is produced using a modified Kaiser window function. 